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A detailed study of a 'PT-symmetric zero-dimensional quartic theory is presented 
and a comparison between the properties of this theory and those of a conventional 
quartic theory is given. It is shown that the PT-symmetric quartic theory evades 
the consequences of the Mermin-Wagner-Coleman theorem regarding the absence 
of symmetry breaking in d < 2 dimensions. Furthermore, the PT-symmetric the¬ 
ory does not satisfy the usual Bogoliubov limit for the construction of the Green’s 
functions because one obtains different results for the h ^ 0~ and the h ^ 0^ limits. 
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INTRODUCTION 


This paper examines the structure of PT-symmetric [T] —93^ theory in d = 0 dimensions 
and makes detailed comparisons between this theory and conventional theory. This work 
is a natural sequel of our studies of PT-symmetric i(p^ theory [2113]. We will see that the 
structure of the theory is much richer and more elaborate than that of conventional 
theory and also of the cubic theories. 

This paper is organized as follows. In Sec. we review the main features of the conven¬ 
tional theory in d = 0 dimensions in the presence of an external linear source h. The 
partition function for this theory is Jdipe~^^‘^\ where V{(p) = h(p -\- /2 -|- /2A. We 

analyze this theory using the method of steepest descents. There is a critical value of mfi\ 
For m? > ml all three saddle points are real, while for one saddle point is real and 

the other two form a complex-conjugate pair. In both cases semiclassical approximations to 
the partition function and the corresponding Green’s functions are obtained. Because d < 2, 
the one-point Green’s function vanishes even when (When d > 2, if 

symmetry breaking occurs.) 


In Sec. Ill the corresponding analysis is performed for the PT-symmetric theory in which 


g and h are replaced by —g and ih. Again, we find a critical value of m^. For > m^, 
one of the three saddle points is purely imaginary and the other two are PT-conjugate. For 
all three saddle points lie on the imaginary-:^ axis. When > m^, new and 
surprising features appear. We hnd that for h > 0 only one saddle point contributes to Z[h], 
while for h < 0 all three saddle points contribute. As a consequence, the partition function 
Z[h] for h > 0 differs from the partition function for h < 0 and the Green’s functions 


*Electronic address: cmb@wustl.edu 
^Electronic address: branchina@ct.infn.it 
^Electronic address: emanuele.messina@ct.infn.it 




2 


obtained by taking the limit h —?• 0“ differ from those obtained in the limit h —)■ O’*". 
Thns, the PT-symmetric qnartic theory evades the properties of the Bogolinbov limit for 
the construction of the Green’s functions. (The Bogolinbov calculation consists of taking 
derivatives of the vacuum persistence amplitude with respect to the external source h, and 
then performing the h —)■ 0 limit, irrespective of the sign of h.) 

When m? < the Green’s functions obtained in the h —)■ O’*" and h —)■ 0“ limits are 
the same. However, unlike the conventional qnartic theory, the one-point Green’s function 
Gi is nonvanishing even though d = 0; the Green’s function has a purely imaginary value. 
The nonvanishing of Gi contradicts the Mermin-Wagner-Goleman (MWG) theorem on the 
absence of symmetry breaking in d < 2 dimensions |1]. We also find that when m? > m'^, the 
one-point Green’s function is nonvanishing. Section |IV| gives a summary and conclusions. 


II. CONVENTIONAL THEORY 

In this section we review the ordinary ip‘^ theory in d = 0 dimensions so that we can 
compare it later with the corresponding PT-symmetric theory. We begin by considering the 
generating functional for the Green’s functions, which is given by the (normalized) integral 

■ + 00 

( 1 ) 

-OO 

The potential is given by 

V{(p) = h(p+^(p'^ + (2) 

where g is assumed to be positive to guarantee the convergence of ([^, but both and h 
can take positive and negative values. The normalization coefficient ^|m^|/(27r) in ([^ is 
chosen so that for the noninteracting theory 



Zfree [O] 



1 . 


There is an implicit factor of 1/h in the exponent in the integral ([^. In this paper we 
calculate the partition-function integral in the semiclassical approximation for which h is 
small. To do so we use the method of steepest descents |5j. We then identify the Green’s 
functions as coefficients of powers of h in the expansion of the partition function Z[h]. 

In order to evaluate the integral in ([^ asymptotically in the limit h —)■ 0, we consider the 
equation that determines the location of the saddle points: 


V'{ip) = h + m^(p + = 0. (3) 

0 

When solving (^, we note that two different situations arise depending on the value of m^. 
Specifically, by defining the critical value = — {9gh‘^) ' /2, we see that the qualitative 
behavior of the three solutions of (|^ depends on being larger or smaller than We 

begin by considering the case 
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A. The m? > rti^ case 


For m? > rril there is one real solution ip^ and there are two complex-conjugate solutions 
and <p_ = to (|^. For small h these solutions are 






9 


^9 


8 V 2mio 


+ 0{h^ 


For h > 0 (h < 0) the solutions form a triangle with one vertex on the negative (positive) 
side of the real-<p axis and the other two vertices in the half-plane Re (p > 0 (Rev? < 0). 

The integral in ([^ is approximated by using Laplace’s method in the limit h —)■ 0. On 
the real axis, there is only one Laplace point [zero of R'(v?)], namely v? 0 ) which is a minimum 
of the potential V{(p). Thus, for all h the partition function Z[h] is approximated as 


Z[h] = exp 




( 4 ) 


Inserting v?o into (|^ and expanding in powers of h, we get the formal Taylor series 


Z[h] = l+{ — 




2m® 




( 5 ) 


As a consequence of the normalization [see ([^] the coefficients of powers of h in (|^ give 
the Green’s functions of the theory. The hrst two connected Green’s functions G\ = Gi and 
GI = G 2 - Gl are 

g; = o, 




2m® 


These expressions are the usual perturbative result (in powers of the coupling constant) [6]. 


B. The < mi case 


For m? = —p? < ml all solutions to (|^ are real: 


l^l — 27r\ /16'|-|-27r 0(±h) \ 

(Po = pcos( -^- 1, (p±=pcosl -^- 1, 


where 0(h) is the Heaviside step function 


1 (h>0), 
0(h) = <; I (h = 0), 
0 (h<0). 


( 6 ) 


and 


/ 2/i^ / / 8/i® 

p = -2sgn(h)^/—, e = -sgn(h) arctan \J-1 + 
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Note that (^_ always lies on the negative-99 axis while 99+ always lies on the positive-99 axis. 
For small values of h in (|^ we hnd that 


^0 


— + 0 {h?), (p± 





+ o{e). 


By examining the second derivative of the potential V"{^p) = m? -\- (799^/2 at these three 
saddle points, we hnd that on the real-99 axis, 99^ are two local minima and that 999 is a 
local maximum of V. For h > 0 the absolute minimum is at 99_ while for h < 0 the absolute 
minimum is at 99_|_. Consequently, the contribution to the asymptotic behavior of Z[h] in 
Q from 99± is given by 


^[OT 


exp 




^ln(V""(99±)//i^) 


Observe that when h > 0 (h < 0 ) the asymptotic behavior of Z[h] is dominated by (p- (99+) 
and that the contribution from 99+ (99-) is subdominant (exponentially small). Expanding 
the above formula in powers of h, we obtain 


m 

Z[ 0 ] 


(1 + ^h,o) 


1 + 



0(h) 



+ e(ft)t 



e(-A)l 



h 


V /x 2 g 32 


y + 0{h^) 


( 7 ) 


This expression is discontinuous at h = 0 . A doubling occurs at h = 0 because at this 
special point there are no subdominant contributions; that is, 99^ contribute equally to the 
semiclassical result ([^ when h = 0. 

Finally, we identify the connected Green’s functions G\ and G2 from ([^: 


G5 = 0, 


1 6/x^ 13^ 

/i2 g 32 /r 6 - 


Note that for d = 0 dimensions, the one-point Green’s function G\ vanishes even when 
As is well known, when d < 2 there is no spontaneous symmetry breaking 
because there is tunneling between the two minima 99_|_ and 99_ of the potential [ 7 ]. In our 
case (d = 0) the vanishing of G^ is due to the cancellation between the two contributions 
to Z[h] coming from the vicinity of 99+ and ip-. Moreover, G2 does not have the usual 
perturbative behavior G2 = ^ + 0 {g). This is because there are contributions from two 
saddle points, which give a typical nonperturbative result [8] . 

Let us compare the above result with the hypothetical d > 2 case for which the volume 
factor suppresses the tunneling (and so there is symmetry breaking). Under this supposition 
we would consider contributions from only one of the two vacua, 99_|_ or ip_. If we were to 


select 99_, then would be replaced by 


z(-)[ 

z(-)[ 


Thus, from ( 7 ) we would hnd that 


i-)Gl 



1 5 g 



( 8 ) 

( 9 ) 
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In (|^ we observe the expected nonvanishing vacnnm expectation value for (p, which is given 
by the tree-level result plus a perturbative correction. Correspondingly, in 

0 we see that the propagator for the fluctuation around has the typical perturbative 
form G^ = ^ + 0{g). 

In the next section we will study the PT-symmetric version of the theory and calculate the 
corresponding Green’s functions. We will discover some intriguing new properties, different 
from those of the conventional theory. In particular, we hnd that even though the formulas 
for the Green’s functions in (|^ and ([^ are only hypothetical, the correct formulas for the 
Green’s functions G1 and G 2 of the corresponding "PT-symmetric theory can be obtained 
by analytic continuation of the formulas in (|^ and ([^. 


III. VT- SYMMETRIC THEORY 


In this section we study the behavior of the massive theory in zero dimensions and 
with negative coupling —g {g > 0). We examine this theory in the presence of a linear source 
term ihip in the potential V{p>)\ 


V{ip) =thip+—ip - —ip . 
To this end, we calculate the partition function 

■ 


( 10 ) 


( 11 ) 


where G is a path of integration in the complex-(y9 plane for which the integral converges. 
In the present case this means that G terminates inside of two Stokes wedges of angular 
opening vr/d, one centered about —vr/d and the other centered about —Svr/d [T]. We also 
calculate the hrst two connected Green’s functions Gl and G^. 

As in Sec. we evaluate the integral in ( |II| ) in the steepest-descent approximation, 
which consists of localizing the integrand around a hnite number of saddle points for which 
the real part of —V is maximal. On the steepest-descent contours ImV is constant. To 
locate the saddle points we must solve the equation 


V\ip) = ih + m^ip — -ip^ = 0. (12) 

6 

Then, we must hnd the two steepest-descent paths in the complex-(y9 plane on which the 
saddle points are the maximum values of —Re V. These paths are also constant-phase 
contours (ImR = constant). It is crucial to identify the constant-phase contour that can be 
deformed into the original integration path G. 





A. The m? > rn^ case 
1. Saddle points 
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Equation (12) has three solutions: 

ifo = -ijd, 


V?+ = + 


(13) 

(14) 

(15) 


where a and (3 are real numbers given by a = 2m?/A + A/g and (3 = 2m? jA — A/g, and 

A = ■\/—2gm'^ sgn(B)jBj^^^ and B =—sgn(h) + \/T — m^/m^. 

Since the theory is VT symmetric, the solution which has a positive-real component, 
is transformed into the (p- solution, which has a negative-real part, by a reflection through 
the imaginary axis. The solution ipo lies on the imaginary axis. We expand (13)-(15) for 
small h to obtain much simpler expressions for the saddle points: 


Po = -i\ + 0{h^), ip± = ±<^ 


16m‘^ h 

-h i 

9 


± — 


2m^ 8 V 


^9 


+ 0{h^ 


Note that for h > 0, and lie in the Im<y9 > 0 half-plane and for h < 0 they lie in the 
Im 9 ? < 0 half-plane; for positive values of h the solution pq lies on the negative-imaginary 
axis and for negative values of h it lies on the positive-imaginary axis. For h = 0 the saddle 
points all lie on the real axis. 

To identify the steepest-descent curves, we let (p = u -t- iv so that the real part of the 
potential in (10) is 


Re V{u, v) = ~ ~ ~ 6u^v^ + u^) 


hv. 


(16) 


Then, from the saddle-point equation (12) we have 


^ T) j/ 2 9 3 I 9 2 n 

—Re V = m u - u + -uv = 0, 

ou 6 2 

Next, we construct the Hessian matrix 


2,9 2 9 3 u n 

—Re V = —m V 3 —u v - v — h = 0. 

ov 2 6 


92 


ReR 


92 


ReR 


U — ( dudu dudv 

^ReR 


dvdu 


dvdv 


m^ — ^{v? — u^) guv 

guv —m? -I- — u^) 


and evaluate it at each of the saddle points (13)-(15). The eigenvectors associated with 
the positive eigenvalues of this matrix provide the steepest-descent directions of —V{u,v). 
On a steepest-descent path the contribution is localized at the saddle point, so we can use 
Laplace’s method to evaluate the integral. For the saddle point po, Uq = 0 and Vq = —(3 [see 
(13)] and thus 


H{p,) = 


m^ + 


0 

-m^ - |uo 
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This matrix is already diagonal. One can easily see that for the eigenvalne 

Aq = + 1^0 is positive, so the u axis is the direction of steepest descent from (po- 

For the other two saddle points and v± = and thus 




-| m | gu±v± \ 
gu±v± ) ■ 


For n = M+ and n = the positive eigenvalue is 

= -ga\/ + 3/3^ (17) 

and the corresponding unit eigenvector is 


wv = 


-,sgn(/?)- 


a/ 2 (4a^ + 3/3^)^/'^ 


a/ 2 (4a^ + 


The negative eigenvalue is 

A/ = —^ga\/Aa^ + 3/^ (18) 

with unit eigenvector 

1 ^ 2 a + ^^2 + 3^2 ^ ^- 2 a + ^^2 + 3/2 

71 (4a2 + 3/2)i/4 (4a2 + 3/2)i/4 






For M = M_ and n = n_ the eigenvalues Aj^ and A 2 and the corresponding unit eigenvectors 

An 


i and W 2 are obtained from the previous ones by making the replacement a —)■ 


-a. 


inspection of (17) and (18) shows that for this case the positive eigenvalue is A 7 = — A^ 
while the negative one is A 2 = —A/. For notational simplicity, in the following we label 
the two components of the unit vector using the letters a and b: wf = {a,b). Also, we 
replace A^ (A 7 ) with A. With this choice of parameters a and b we also get tiJ/ = (— 6 , a), 
Wi = (—a, b), and W2 = {b, a). As explained previously, for our calculation we are interested 
in the positive eigenvalues and the corresponding eigenvectors. With the notation we have 
introduced the positive eigenvalue is A in both cases {ip+ and ip_) and the corresponding 
eigenvectors are (a, 6 ) and {—a,b). 


2. Constant-phase contours 


We now determine the constant-phase contours in the (p plane that are also the steepest- 
descent paths from the saddle points. Our goal is to deform the original path C in ( 11 ) 
into a new set of steepest-descent contours passing through the saddle points (13)-(15). As 
explained earlier, on this new path ImlA is constant (or, more precisely, piecewise constant), 
and the integral in ( [IT| is approximated by a sum of gaussian integrals localized at the 
saddle points. 

The constant-phase contours are determined by the equation 


ImV 


hu + m?uv — 


1 3 1 3 

-gu V -\—guv = c. 
6 6 


(19) 
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FIG. 1: [Color online] Constant-phase contours in the {u,v) plane for /i > 0, with the arbitrary 
choice m? = 0.1, g = 0.2, and h = 5.0. The constant-phase contours associated with the saddle 
point (/9o are green, those associated with (/?+ are blue, and those associated with are red. The 
steepest-descent paths associated with these three saddle points are solid, and the steepest-ascent 
paths emanating from the saddle points are dashed. The arrows indicate the down directions. 
The original contour terminates in Stokes wedges in the southwest and southeast quadrants of the 
complex plane. Therefore, for the case h > 0 the only relevant saddle point is po, and the original 
integration path must be deformed into the contour labeled C. For this case the saddle points ip± 
play no role in the asymptotic evaluation of the partition function. 


The constant c will take the valnes c+ , cq, or c_ when we evalnate (19) at ip+, (po, or 
respectively. We now stndy the two cases h > 0 and h < 0 in tnrn. 

The h> 0 case: We first consider the case h > 0. Fignre shows the saddle points (po, T+y 
and (p_ of ([l^-([l^ in the {u, v) plane for a specific (althongh irrelevant) nnmerical choice 
of the parameters together with the constant-phase contonrs passing throngh these points. 
The solid lines are the steepest-descent (constant-phase) contonrs associated with the three 
saddle points. The phases are c+, c_, and cq, respectively. 

In general, to evalnate the partition fnnction Z[h] of ( [IT| in the saddle-point approxima¬ 
tion, we mnst identify the contonr C that passes throngh one or more of the saddle points 
and that at the same time is a steepest-descent path that terminates inside the Stokes 
wedges in the sonthwest and sontheast qnadrants of the complex plane. An inspection of 
Fig. [^shows that when h > 0, of the three saddle points only one of them, namely cpo, fnlhlls 
these conditions and that the solid cnrve labeled C that passes throngh ipo is the desired 
contonr. For ip = (pQ = wq we have \mV{ipo) = 0 [see (19)], so the constant-phase contonr 







9 


C is given by (19) with c = 0; thus, the equation for the curve C is 

Qh Qm^ 

- 1 - 

gv 9 


+ = 0 . 


This curve terminates in the middle of the Stokes wedges. By evaluating the integral in ( pT| 
along this constant-phase contour C, we obtain 


Z[h] = 


\m^ 


271 


dip e 


'c 


-v(^) ^ r du(l + ] e.-^-v[uMu)] 


277 


du 


Expanding the real part of the potential around the saddle point {uq = 0, vq = —(3) to 
second order in {u — uq) and {v — vo), we get [with Re V (<po) = Re Vq] 

Re V[u, n(M)] ~ Re Rq + ^ (« - ^o, - -^ 0 ) • H{(po) ■ ^ . 


Therefore, the saddle point approximation for Z[h] is 

Z[h] = exp (^-Re Rq - ^ hi , 

where Re Rq = —m^/3^/2 — gl3‘^/24: + h(3 and Aq = 

We expand the partition function Z[h] in (20) in powers of h (/3 ~ and get 


( 20 ) 


Z[h] = 1-1 —+ 




2m® 


-^ + o{h^). 


( 21 ) 


From this result we can obtain the one-point and two-point Green’s functions Gi = 


; 1 dZ[h] 
Z[0] dh 


h=0 


d'^Z[h] 

and G2 [i) ^^2 


h=0 


. Finally, the connected Green’s functions are 


= 0 , 

Gi = — + 


2m® 


( 22 ) 

(23) 


We emphasize that despite the negative sign of the quartic coupling constant (—< 0) 
(which one might interpret as implying instability due to unboundedness below), the function 
G 2 is real and positive as a consequence of VT symmetry. It is also worth noting that in this 


case the partition function (20) and the two-point Green’s function (23) has the form one 


would expect from a perturbative expansion in powers of the coupling constant g. In other 
words, the saddle-point approximation to Z[h] gives, as in the conventional positive-coupling 
constant theory, the usual perturbative expansion. For this reason, we call the saddle point 
ipo the perturbative saddle point. 


Let us compare the partition function (21) of the "PT-symmetric theory with the corre¬ 
sponding partition function ([^ of the conventional theory for the m^ > m^ case. Note hrst 
that the potential (10) of the PT-symmetric theory is obtained from the potential ([^ of the 
conventional theory by making the replacements g —)■ —g and h —)■ ih. By making the same 
replacements in the partition function ([^ of the ordinary theory, we obtain the partition 
function (21) of the PT-symmetric theory. Thus, by analytically continuing the partition 
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FIG. 2: [Color online] Constant-phase contours in the (u, v) plane for the case h < 0 with the 
arbitrary choice = 0.1, g = 0.2, and h = —5.0. As in Fig. the constant-phase contours 
associated with the saddle points (po, (p+, and (p- are colored green, blue, and red. The steepest- 
descent paths associated with these three saddle points are solid, and the steepest-ascent paths 
emanating from the saddle points are dashed and the arrows indicate the down directions. As 
in the case h > 0, the original integration contour terminates in Stokes wedges in the southwest 
and southeast quadrants of the complex plane. However, unlike Fig. all three saddle points 
are relevant. The original integration path must be deformed into three constant-phase contours 
labeled C-, Cq, and (7+, which are joined end-to-end. 


function of the conventional theory to imaginary values of h and negative values of g, the 
partition function of the corresponding "PT-symmetric theory is obtained. We will see in 
the following that this simple connection between the two theories is lost when the partition 
function of one of the two theories gets contributions from more than one saddle point while 
the partition function of the other theory gets a contribution from only one saddle point. 

The h < 0 case: When h < 0, the saddle points (po, T+i (p_ and the constant-phase 
curves are obtained from those of the h > 0 case by reflecting about the u = 0 axis (compare 
Figs. and [^. This has crucial consequences for the saddle-point approximation of Z[h] in 
(0. Indeed, Fig. shows that the contour needed to evaluate the integral is quite different 
from the contour C of Fig. [T} As shown in Fig. the integration path is now obtained by 
joining three constant-phase contours, C-, Co, and C+. 

For h < 0 (unlike the h > 0 case) all three saddle points contribute to the integral in 
p. Thus, along the path C_ -|- Co -|- C+ we decompose the integral as follows: 

Z|;i] = yiP (/ + j + J (l/pe-'''-*'' 



(24) 
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We consider each of these integrals in turn: For the integral along Co, the calculation is the 
same as for the h > 0 case, and we obtain 




d(p e 




exp 


'Co 


-ReW 






(25) 


The curve C+ is given by (19), and using (14) we get 


a/s ^ a/s 2 n 3 


9 


c = c+ = Iml/((p+) = h—a + —m a/? - 




tie 


We parametrize C+ in terms of the variable u = u^{t) and v = n+(t) and note tha 
vector wf = {a,b) is the unit vector tangent to C+ in (p+ = u+ + iv+ = ^a + i^l3 [see (14) 
Expanding u+(t) and n+(f) to hrst order in t, we get u+(t) — ~ at, v+(t) — ~ bt, so 

the expansion of Re V{ip) around (p+ to second order in t is 


ReR(^) = ReR((p+) + - 


[u+{t) - u+,v+{t) - v+] ■ H{p>+) ■ 


— ReR((p+) + 

where we have used + 6^ = 1. The integral along C+ is then approximated by 

dp> ~ (a + ib) exp 




271 


IC+ 


-Re V((p+) - I In ^-ic+ 
2 


(26) 


Last, we calculate the contribution to Z[h] in (24) coming from the integral along the curve 
C_. Inserting (15) into (16) and (19), we hnd that c_ = —c+ and ReR(</9_) = ReR(</9+). 
Around ip_, the curve is approximated as u_{t) — m_ ~ —at and v_{t) — r>_ ~ bt. 
Following the same steps as for the evaluation of the integral along C+, we hnd that 


|m^| 


d(p e 


-v{v) 


ib) exp 


<C- 


-Rel/((p_) - ^In A 


*c_ 




(27) 


Combining (25), (26), and (27), and using the notation c = c+ = —c_ and ReR = 
ReR((p+) = ReR((p_), the hnal result for the saddle-point approximation to Z[h] is 


Z[h] = exp ( —Re Vq 


bnA 




-|- 2 (a cos c -|- 6 sin c) exp 


-ReR-Jln^ 
2 


(28) 


The contribution to Z[h] for h < 0 that appears in the hrst term of (28) agrees with the 
result for Z[h] in (20) for h > 0. However, in (28) an additional term appears, which comes 
from the contribution of the two additional nontrivial saddle points (p+ and (p_. These 
saddle points do not contribute to Z[h] in the h > 0 case. We emphasize that the h > 0 
and the h < 0 cases are diherent because the saddle-point contour passes only through the 
ifo saddle point when h > 0 (see Fig. [^, while the saddle-point contour passes through all 
three saddle points, (p-, ipo, p>+, when h < 0 (see Fig. |^. Thus, the h > 0 and the h < 0 
cases cannot be obtained from one another by simply changing the sign of h in the partition 
function. 
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The peculiar discontinuity in the partition functions for the h > 0 case (20) and the h < 0 
case (28) only occurs in the VT symmetric theory and does not appear in the conventional 
theory. As a consequence, the Green’s functions obtained from (20), which are the limit 
as h —?• 0^ of derivatives of Z[h], are different from the Green’s functions obtained from 
(28), which appear in the limit h —)■ 0“. To see this difference explicitly, we expand (28) in 
powers of h. Keeping only terms up to 0{h‘^) we get 


Z[h] _ f 1 _ / 1 9 \h^ 

-^[0] \V ^ V 2m^J ^ 2m®/ 2 ’ 

from which we read off the connected Green’s functions and G^'- 


G5 


• f [IL\ 

yy g Sm? V 2m? j 

1 9 

-^- —■ 

m^ 2m® 



(29) 

(30) 


Note that the two-point Green’s functions G 2 in (23) and (30) for the cases h > 0 and h < 0 
are the same. However, the one-point Green’s functions G\ in (22) and (29) are different. 
This difference is a violation of the usual Bogoliubov limit [9] for constructing the Green’s 
functions. For the conventional theory, the Green’s functions are obtained as the h —)■ 0 
limit of Z[K\ derivatives, independently of the sign of h. 


B. The case 

Let us consider now the evaluation of Z[h] in ( |TT| ) for the case m^ = — /i^ < m^ (where 
/i^ > 0). As before, we first look for the saddle points and then determine the constant-phase 
contours for the saddle-point approximation of Z[h]. 


1. Saddle points 


The saddle points are given by the solutions to (|12|) and are purely imaginary: 

'|0| - 27r' 


(p± = ip cos 


|0| +27r0(±h)\ 

3 ) 


Pq = ip cos 


(31) 


where 


' 2 /i 2 


p = —2sgn(h)W-, 6 = —sgn(h) arctan 

9 


-1 + 


8 p® 


For any h, lies on the negative-imaginary axis in the complex-ip plane, (p+ lies on the 
positive-imaginary axis, and depending on the sign of h, (po lies on the positive-imaginary 
or negative-imaginary axis. Expanding the saddle-point solutions in (31) in h, we get 


h 


6 /i^ 


h 




p- 


2 /i 2 


^9 

2/iio 8 


+ 0{h^ 
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Next, letting = u + iv, we find that the real part of the potential (10) is 


Re V{u, v) = 


~{u^ — v^) — ~ + v'^) — hv. 


We then calculate the Hessian matrix, which at the three saddle points (31) is diagonal: 


H = ( 


dvdu 


ReV 

ReV 


duq, 

dvdv 


ReV 

ReV 


-/i 


2 + 1^2 
0 


0 


/i" 


2 y 2 
2 ^ 


(32) 


where n is uq = i(po or v± = (p±, depending on which of the saddle points we consider. 

As before, the local directions of the steepest-descent paths from each of the saddle points 
are determined by hnding the eigenvectors of the Hessian matrix with u and v at the saddle 
points. Because H has a diagonal form, the eigendirections are parallel to the u and the v 
axes. The signs of the two eigenvalues —jV + and jV — are found by substituting 


n+, n_, and Vq [see (31)] in (32). We hnd that —is negative for ipQ and positive for 


and </9_. The opposite is true for pr — 


2. Constant-phase contours 


To hnd the constant-phase contours we must solve the equation 

ImH = hu — jjc^uv — gu^v/6 guv^/6 = c. 


(33) 


where V is the potential in (10) and the constant c is determined by evaluating ImH at 


each of the saddle points ip+, ipo, and ip-. The real part of all three saddle points vanishes. 


so from (33), we have c = 0 for each saddle point. Factoring (33) gives a linear equation. 


namely u = 0, and a cubic equation for the constant-phase contours: 

— [Qfi^/g u^) V + 6h/g = 0. 


(34) 


The u = 0 contour is just the v axis and passes through all three saddle points (p+, (p_, 
and (po- The cubic equation (34) provides the remaining three constant-phase contours, one 


for each of the three saddle points. The constant-phase contours are plotted in Fig. [^for 
a specihc choice of the parameters. In the left (right) panel we consider the h > 0 (h < 0) 
case. Observe that in both cases, the only path that terminates in the Stokes wedges is the 
path that passes through (p_. This path is labeled 0_. 

Laplace’s method along the path 0_ approximates Z[h] as exp Re R_ — | In , where 

ReVT = Rel/((p_) = ia^v^_/2 - gvt/24: - hv_ and A_ = 
normalized partition function expanded in powers of h is 




-|- gv^/2 > 0. Thus, the 


m 

Z[0] 


= 1 




g 


4 V 




9 


32/1^ 




This equation determines the hrst two connected Green’s functions: 




G c _ 

2 “ 


2/i2 


16 yU® 


( 35 ) 
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FIG. 3; Constant-phase contours in the {u,v) plane for the potential V in (10) for the parameters 


m 


2 _ 


= —0.1, g = 0.6. The left panel portrays the case /i > 0 for the specific value h = 0.02 and 
the right panel portrays the case h < 0 for the specific value h = —0.02. One contour lies on 
the imaginary axis and the other three contours emerge horizontally from each of the three saddle 
points (/9+, (p-, and (pQ. The only relevant contour is the one that passes through because this 
contour terminates in the Stokes wedges, which lie in the southwest and southeast quadrants of 
the {u, v) plane. This contour is labeled C_ and is represented as a solid line. All other contours 
are indicated by dashed lines. 


Note that the connected two-point Green’s function G 2 is precisely what one obtains from 
perturbation theory. In particular, the squared mass of the fluctuation above the vacuum 
is = —2m^, as in the conventional theory. In contrast, Gi is nonvanishing. Evidently, 
the PT-symmetric theory evades the MWC theorem |3] (according to which there cannot 
be spontaneous symmetry breaking below d = 2 dimensions). 


IV. SUMMARY AND CONCLUSIONS 

We have studied the zero-dimensional Euclidean PT-symmetric (p^ theory in the pres¬ 
ence of an external source h and have compared it with the corresponding conventional 
theory. We calculated the partition function and the Green’s functions using the saddle- 
point approximation. As is the case for the conventional theory, the "PT-symmetric theory 
possess distinct phases depending on whether the squared mass rn? is greater than or less 
than the critical value These phases are characterized by different Green’s 

functions. Furthermore, when the PT-symmetric theory exhibits two distinct 

subphases depending on whether h —)■ O’*" or h —)■ 0“. 

The subphases can be described as follows: For and for h > 0, the partition 

function and the Green’s functions of the PT-symmetric theory can be obtained from the 
corresponding functions of the conventional theory by analytic continuation to imaginary 
values of the external source h (h —)■ +ih) and to negative values of the coupling g {g ^ ~9)- 
This connection between the two theories arises because in both cases the partition function 
receives contributions from one saddle point only; the saddle point of the PT-symmetric 












15 


theory can be obtained from the saddle point of the conventional theory by an anticlockwise 
rotation of 7 r /2 in the complex-(p plane. On the other hand, for and h < 0, while 

the partition fnnction of the conventional theory is still dominated by the pertnrbative saddle 
point, the partition fnnction of the "PT-symmetric theory receives contribntions from the 
rotated perturbative saddle point and also from two nonperturbative saddle points. This 
subphase behavior of the PT-symmetric theory is a new and unexpected feature. The 
Green’s functions obtained by performing the two limits h —)■ 0+ and h —)■ 0“ are different, 
and this is a clear violation of the usual Bogoliubov calculation of the Green’s functions. 

The PT-symmetric theory differs from the conventional theory in other respects. In 
addition to the behavior of the phase, new phenomena appear in the < ml 

phase. Spontaneous symmetry breaking in the conventional quartic theory is forbidden for 
d <2 due to tunneling between the degenerate minima; this absence of symmetry breaking 
is due to the MWG theorem |3]. The MWG theorem certainly holds for d = 0 and it 
implies that for the conventional theory the one-point Green’s function G\ vanishes even 
when mf < mf. In the context of the saddle-point approximation, the vanishing of 
is due to the cancellation of contributions to the partition function coming from the two 
nonperturbative saddle points. Thus, the PT-symmetric theory evades the MWG theorem 
in much the same way that supersymmetry evades the Goleman-Mandula theorem, namely, 
by circumventing the assumptions needed to prove it. 

However, for the PT-symmetric quartic theory the dominant contribution to Z[h] comes 
from just one saddle point. Thus, the < mf phase of the PT-symmetric theory is 
characterized by a nonvanishing Gi when d = 0. Hence, we would expect that even when 
0 < d < 2, in the PT-symmetric (p‘^ theory a phase transition will occur that is not triggered 
by the usual spontaneous-symmetry-breaking mechanism. Moreover, we have demonstrated 
the surprising result that while Gi = 0 in the conventional theory, Gi 7 ^ 0 when > ml 
in the PT-symmetric theory. 

Let us make a further comparison between the conventional and the PT-symmetric quar¬ 
tic theories. We have seen that when > ml and h < 0, the partition function of the 
PT-symmetric theory is not an analytic continuation of the corresponding partition function 
of the conventional theory. This is because the PT-symmetric Z[h] receives contributions 
from three saddle points, while only one saddle point contributes to the partition function 
of the conventional theory. In the < ml case, the situation is similar although the roles 
of the ordinary and the PT-symmetric theories are exchanged. The partition function of 
the conventional theory receives contributions from two saddle points, while Z[h] of the 
corresponding PT-symmetric theory gets contributions from only the one saddle point (p_ 
(see Fig. [^. We stress that a cancellation leading to Gi = 0 occurs when there are two 
saddle points, but this is impossible when there is only one saddle point. 

We have pointed out in Sec. in our discussion of the conventional theory that we can 
artihcially induce spontaneous symmetry breaking by discarding the contribution to Z[h] 
from ip^. For d <2 this procedure is artihcial because both the left and right saddle points 
contribute to Z[h] so the theory cannot exhibit spontaneous symmetry breaking. Of course, 
one of the two saddle points must be discarded only when d > 2. We have entertained this 
possibility for the d = 0 conventional theory in Sec. [TTJonly to compare the Green’s functions 
in (35) with the analytic continuation (a rotation in the complex-(y9 plane) of the Green’s 
functions in ([^ and (|^. 

To conclude, the structure of the quartic PT-symmetric theory in zero dimensions is 
much richer than that of conventional quartic theory. It is also far more elaborate than that 
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of the cubic theory EE]. Our immediate objective now is to extend the results for the 
■PT-symmetric quartic theory to the d > 0 case, just as we did for the cubic theory. 
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